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ABSTRACT 

We report the spectroscopic confirmation of four further white dwarf members of Prae- 
sepe. This brings the total number of confirmed white dwarf members to eleven making this 
the second largest collection of these objects in an open cluster identified to date. This number 
is consistent with the high mass end of the initial mass function of Praesepe being Salpeter in 
form. Furthermore, it suggests that the bulk of Praesepe white dwarfs did not gain a substan- 
tial recoil kick velocity from possible asymmetries in their loss of mass during the asymptotic 
giant branch phase of evolution. By comparing our estimates of the effective temperatures 
and the surface gravities of WD0833H-194, WD0840H-190, WD0840H-205 and WD0843H-184 
to modern theoretical evolutionary tracks we have derived their masses to be in the range 
0.72— O.76M0 and their cooling ages ^300Myrs. For an assumed cluster age of 625±50Myrs 
the infered progenitor masses are between 3.3— 3.5M0. Examining these new data in the con- 
text of the initial mass-final mass relation we find that it can be adequately represented by 
a linear function (ao=0.289±0.051, ai=0.133±0.015) over the initial mass range 2.7M0 to 
6M0 . Assuming an extrapolation of this relation to larger initial masses is valid and adopting a 
maximum white dwarf mass of 1 .3M0, our results support a minimum mass for core-collapse 
supernovae progenitors in the range ~6.8-8.6M0. 
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1 INTRODUCTION 

The initial mass-final mass relation (IFMR) characterises the 
amount of material stars with primordial masses M^1OM0 cast 
out into interstellar space during post main sequence evolution, en 
route to becoming white dwarfs. Accordingly, the form of this rela- 
tion is of considerable importance to investigations relating to the 
chemical evolution of the Milky Way and galaxies in general. The 
details of the upper end of the IFMR also have relevance to su- 
pernovae studies since theoretical predictions of the rate of these 
explosions are rather sensitive to the assumed minimum mass of a 
core-collapse progenitor. Furthermore, the form of the IFMR con- 
veys information on the mass loss processes which occur during 
the final stages of stellar evolution, which are difficult to model in 
a physical context. 

Arguably the most robust method with which to place (semi-) 
empirical constraints on the form of the IFMR is via the study of 
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white dwarfs in open clusters (e.g. Weidemann 1977, Romanischin 
& Angel 1980, Weidemann 2000). Since the constituents of an open 
cluster have a common age, determinable from the main sequence 
turn-off mass (e.g. King & Schuler 2005), the lifetime of a pro- 
genitor star can be estimated from the difference between this age 
and the cooling time of the resulting white dwarf. Subsequently, the 
progenitor masses can be determined by comparing their estimated 
lifetimes to the predictions of stellar evolutionary models. 

Unfortunately, until quite recently, the small numbers of WDs 
recovered in each open cluster (~30 WDs in 14 clusters), their in- 
trinsic faintness and the significant distances involved has meant 
that the uncertainties in clusters ages and in white dwarf mass de- 
terminations (and hence cooling time estimates), have conspired to 
produce significant scatter in the semi-empirical IFMR (e.g. Claver 
et al. 2001). Encouragingly, with time on 8/lOm class telescopes 
now readily available, this situation is beginning to improve. For 
example, Kalirai et al. (2005) and Williams et al. (2005) have spec- 
troscopically identified 16 and 6 likely white dwarf members of the 
rich but relatively distant open clusters M37 and M35 respectively. 
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Figure 1. The results of our fitting of synthetic profiles (thin black lines) to the observed Balmer lines, H-/3 to H-8, (thick grey lines). The flux^ units are 
arbitrary. 



These two populations of white dwarfs are found to be consistent 
with a monotonically increasing IFMR for progenitor mass ranges 
M~2.8-3.4Mo and M~4.5-5.5M0, respectively. 

Nevertheless, while significant additional pieces of the jigsaw 
puzzle have recently been put in place, much work remains to be 
done before we are able to claim a thorough understanding of the 
IFMR. Accordingly, wc recently embarked on a search for white 
dwarfs members of the modestly rich and relatively nearby Prae- 
sepe open cluster (M44, NGC2632). Using proper motions mea- 
surements and photographic photometry drawn from the USNO- 
Bl.O catalogue and SuperCosmos archive we identified 6 new can- 
didates in a 5° x 5° region centered on the cluster, to add to the 
five previously known Praesepe white dwarfs. Futhermore, two of 
these new candidates were confirmed spectroscopically. Full details 
of these objects are given in Table 1 of Dobbie et al. (2004). 

In this work we present new low resolution optical spec- 
troscopy which confirms that the four remaining candidates are also 
white dwarfs. Additionally, we use new data to provide refined esti- 
mates of the effective temperatures and surface gravities of LB5959 
and WD0837+218. We re-examine the membership status of each 
of these objects and present a new astrometric measurement of the 
high mass "Praesepe" white dwarf LB5893 which establishes that 
it has a proper motion which is consistent with cluster membership. 
We discuss the implications of our findings in relation to the initial 
mass function of Praesepe and kinematic effects of mass loss in the 
final stages of the evolution of the white dwarf progenitors. Adopt- 
ing an age of 625±50Myrs for the cluster, we estimate the masses 
of these progenitor stars and examine our results in the context of 
previous work on the initial mass-final mass relation. 



OPTICAL SPECTROSCOPY OF THE REMAINING 
FOUR CANDIDATE WHITE DWARF MEMBERS 



We have acquired optical spectra (3200-6000A) of the four re- 
maining candidate white dwarf members of Praesepe from Dobbie 
et al. (2004), using the WiUiam Herschel Telescope and the ISIS 
spectrograph. The observations were conducted during the service 
time nights of 2005/01/10 and 2005/1 1/07 and the visitor night of 
2006/02/02. Sky conditions were good on the first and third nights 
with clear skies and seeing ~0.7-0.9". On the second night there 
was some thin and patchy cirrus and seeing was slightly poorer, 
~1.0-1.2". All data were obtained on the blue arm of ISIS, using 
the EEV12 detector, the R300B grating and a slit width of 1" to 
provide a spectral resolution of « 3.5A. The total exposure time 
for each object was 45 minutes made up from a series of 15 minute 
integrations. The CCD frames were debiased and flat fielded using 
the IRAF procedure CCDPROC. Cosmic ray hits were removed us- 
ing the routine LACOS SPEC (van Dokkum 2001). Subsequently 
the spectra were extracted using the APEXTRACT package and 
wavelength calibrated by comparison with the CuAr-i-CuNe arc 
spectra. To facilitate the removal of the instrument signature from 
the science spectra we obtained observations of the spectral stan- 
dard stars PG0843+546 (Massey et al. 1988) and G191-B2B (Oke 
1990). 
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Table 1. Details of the four new confirmed white dwarf candidate cluster members (top) 
and the two white dwarfs members identified by Dobbie et al. (2004; bottom) for which 
we have obtained improved spectroscopic data. Masses and cooling times for each star 
have been estimated using the mixed CO core composition "thick H-layer" evolutionary 
calculations of the Montreal Group (e.g. Fontaine, Brassard & Bergeron 2001). 



WD IDinD04 T<,fl(K) log g M(Mq) Tc(Myrs) 

8.18l:JJ;°3 0.72 ±0.02 262^14 

0840+190 candidate 4 14765^1^^ 8.2ll° °^ 0.74 ±0.02 2881^3 

•24j;o;°4 0.76 ±0.03 316t^g 

0843+184 candidate 6 14498l2og S-'^'^toM 0.75 ±0.03 3081^9 



•241^;;;^ 0.76 ±0.04 303125 
0837+218 candidate 3 I6833I272 S.SQIqqj 0-86 ± 0.02 267^^3 



D04: Dobbie et al. (2004) 

3 ANALYSIS OF THE DATA 3.2 Determination of effective temperatures 



3.1 Model wliite dwarf spectra 

The broad hydrogen B aimer lines evident in the data presented 
in Figure 1 are consistent with all four objects without prior 
spectroscopic data being DA white dwarfs. We have used the 
latest versions of the plane-parallel, hydrostatic, non-local ther- 
modynamic equilibrium (non-LTE) atmosphere and spectral syn- 
thesis codes TLUSTY (v200; Hubeny 1988, Hubeny & Lanz 
1995) and SYNSPEC (v48; Hubeny, I. and Lanz, T. 2001, 
http://nova.astro.umd.edu/ 1 to generate a new grid of pure-H syn- 
thetic spectra covering the Tea and surface gravity ranges 10000- 
34000K and log g=7.0-9.0 respectively. We have employed a model 
H atom incorporating the 8 lowest energy levels and one superlevel 
extending from n=9 to n=80, where the dissolution of the high ly- 
ing levels was treated by means of the occupation probability for- 
malism of Hummer & Mihalas (1988), generalised to the non-LTE 
situation by Hubeny, Hummer & Lanz (1994). All calculations in- 
cluded the bound-free and free-free opacities of the H~ ion and 
incorporated a full treatment for the blanketing effects of HI lines 
and the Lyman —a, —(3 and —7 satellite opacities as computed by 
N. Allard (e.g. Allard et al. 2004). In contrast to the grid of mod- 
els used in our previous work where radiative equilibrium was as- 
sumed (Dobbie et al. 2004), these latest calculations include, where 
appropriate, a treatment for convective energy transport according 
to the ML2 prescription of Bergeron et al. (1992), adopting a mix- 
ing length paramater, a=0.6. During the calculation of the model 
structure the hydrogen line broadening was addressed in the fol- 
lowing manner: the broadening by heavy perturbers (protons and 
hydrogen atoms) and electrons was treated using AUard's data (in- 
cluding the quasi-molecular opacity) and an approximate Stark pro- 
file (Hubeny, Hummer & Lanz 1994) respectively. In the spectral 
synthesis step detailed profiles for the Balmer lines were calculated 
from the Stark broadening tables of Lemke (1997). 



and surface 
gravities 

As discussed in Dobbie et al. (2004), comparison between the mod- 
els and the data is undertaken using the spectral fitting program 
XSPEC (Shafer et al. 1991). XSPEC works by folding a model 
through the instrument response before comparing the result to the 
data by means of a x^^statistic. The best fit model representation 
of the data is found by incrementing free grid parameters in small 
steps, linearly interpolating between points in the grid, until the 
value of is minimised. Errors in the Tcft S and log g s are calcu- 
lated by stepping the parameter in question away from its optimum 
value and redetermining minimum until the difference between 
this and the true minimum corresponds to Icr for a given number 
of free model parameters (e.g. Lampton et al. 1976). 

Given the probable age of the Praesepe cluster (~600- 
700Myrs), a number of the white dwarf members may have ef- 
fective temperatures approaching that at which the H-Balmer lines 
reach their maximum equivalent width (Tcft '^12500-13500K). Pre- 
vious spectroscopic studies of DA white dwarfs in this temperature 
regime indicate that the sensitivity to surface gravity of the equiv- 
alent widths of the lower order Balmer lines (e.g. H-/3, H-7) is re- 
duced here (e.g. Daou et al. 1990). Therefore, in the present analy- 
sis all lines from H-/9 to H-8 are included in the fitting process. Fur- 
thermore, given that we assign each line an independent normalisa- 
tion parameter, there is the potential for obtaining two solutions for 
Tcff in this regime (e.g. Gianninas et al. 2005). Hence minimum 
has been approached from both the high and low temperature ends 
of the model grid. 

No convincing fit with Tcff^l3500K is found to any of the 
datasets; the overall spectral shape of the data (ie. lines and contin- 
uum) are well matched by models at the temperature and surface 
gravity solutions determined here. The results of our fitting proce- 
dure are given in Table 1 and shown overplotted on the data in Fig- 
ure 1. As a check for potential systematic errors (e.g. Napiwotzki 
et al. 1999) we have validated the results with an independent anal- 
ysis technique (FITSB2; Napiwotzki et al. 2004) and a different 
set of model atmospheres (the LTE atmospheres used by Koester 
et al. 2001). Differences between individual results were found to 
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Figure 2. A schematic plot of the Praesepe cluster showing stars down to V 
~ 9 and the areas surveyed by Anthony-Twarog (1982, 1984; solid outline) 
and Claver et al. (2001; grey shading). The region included our investiga- 
lion is outUned (dashed grey Une). All objects listed in Table 1 of Dobbie et 
al. (2004; open circles) and the five 'original' white dwarf cluster members 
(open stars) are also overplotted. The locations of the four new spectro- 
scopically confirmed white dwarf candidate members are highlighted (open 
triangles). 



be well within the error limits quoted in Table 1 and no significant 
systematic discrepancies were apparent. Nevertheless, it should be 
noted that the parameter errors quoted here are formal Icr fit errors 
and may underestimate the true imcertainties. 



4 DISCUSSION 

4.1 White dwarf membership status, the initial mass function 
of Praesepe and recoil kicks 

We have used modem evolutionary tracks supplied by the Mon- 
treal group (e.g. Fontaine, Brassard & Bergeron 2001) to determine 
the masses and cooling times of our four new white dwarfs. We 
have adopted the calculations which include a mixed CO core and 
thick H surface layer structure, which make this work consistent 
with other recent studies in this area (e.g. Liebert et al. 2005a). The 
masses and cooling times shown in Table 1 have been derived us- 
ing cubic splines to interpolate between points in this grid. At the 
present level of precision these mass determinations are not sensi- 
tive to our choice of core composition. However, if instead we had 
adopted thin H-layer models these estimates would be systemati- 
cally lower by 0.02Mq . Nevertheless, we can conclude that all four 
white dwarfs have masses which are significantly larger than those 
typical of field white dwarfs, the distribution of which is foimd to 
be strongly peaked at O.565M0 (e.g. Liebert et al. 2005b, Berg- 
eron, Liebert & Fulbright 1995, Marsh et al. 1997). These com- 
paratively high masses and the projected spatial distribution of the 
white dwarfs (Figure 2) argue that our objects are associated with 
the Praesepe open cluster. 



Table 2. Progenitor lifetimes and corresponding masses based on the 
white dwarf cooling times shown in Table 1, the Z=0.019 stellar evo- 
lutionary models of Girardi et al. (2000) and an assumed cluster age of 
625±50Myrs. 
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In an effort to quantify the level of contamination by field 
white dwarfs in our assembly of probable Praesepe members, we 
have repeated our original survey procedure on four 5° x5° fields 
flanking the cluster to the NE, NW, SE and SW. Only 2 objects 
in total have been flagged as candidate white dwarf members of 
Praesepe according to the selection criteria we applied in Dobbie 
et al. (2004). There is no guarantee that spectroscopic data would 
confirm either of these two objects to be a white dwarf. Thus we 
conclude that contamination by field white dwarfs is at a very low 
level. 

Nevertheless, it is worth recalling here that LB5893, one of 
the five "original" white dwarf members of Praesepe, was not re- 
covered by our survey. Since Luyten (1966) measured a proper 
motion of jiaCOS (5=-34 mas yr^^, /X5=-14 mas yr~^, albeit with 
large uncertainties, in Dobbie et al. (2004) we concluded that the 
USNO-Bl.O astrometry of /iacos (5=— 56 mas yr~^, /X5=— 14 mas 
yr~^, may have been adversely affected by the proximity of the star 
KW195. Claver et al. (2001) previously reached a similar conclu- 
sion on finding that their proper motion measurement of this object, 
which relied on POSS-I plates for epoch 1 , also failed to confirm 
cluster membership. Given the lingering uncertainties in the proper 
motion of LB5893 and the rather peculiar location of this object 
in the semi-empirical IFMR (see Figure 11 of Claver et al. 2001), 
we have used the POSS-II F band image (1989/11/08) and a V 
band image (2001/02/04) obtained from the Canada-France-Hawaii 
Telescope archive at the Canadian Astrophysical Data Center to ob- 
tain a refined astrometric measurement for LB5893. Our estimate 
of /IqCos 5=— 34±9 mas yr~^, /i5=— 18±12 mas yr~^ supports 
the original proper motion determination of Luyten (1966) and ar- 
gues strongly that LB5893 is a member of Praesepe. Therefore in 
the subsequent discussion we accept that the number of Praesepe 
white dwarfs is at least eleven. 

As a consequence of the optical (O-E^O) and near-infrared se- 
lection criteria (no 2MASS detection or blue) our survey is biased 
against the detection of white dwarfs in unresolved binary systems 
with stars with masses M^O. IM© (spectral types earlier than M6- 
7). Thus in the context of the cluster simulations of Williams et 
al. (2004), in reflection, it seems more appropriate to compare our 
results to the predicted observable number of single white dwarfs 
(ie. the predicted total observed number - the predicted number ob- 
served in binaries). We find that the observed number is inconsis- 
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Figure 3. Initial mass-final mass for the white dwarf members of the Hyades (open triangles), Praesepe (black circles), M35 (open diamonds), NGC2516 
(open '+'s) and the Pleiades (open star). Sirius B is also shown (open circle). A linear fit to the data based on CO (soUd line)/C(dashed line) cores and the 
relation of Weidemann (2000; dotted line) are overplotted. 



tent with these predictions, for any reasonable value of the maxi- 
mum mass of a white dwarf progenitor (Merit), if a steep power- 
law (r=2) shape is assumed to describe the initial mass function 
(IMF) of Praesepe. In constrast, we find that the observed number 
is consistent with the results of the simulations if instead the IMF 
is assumed to have been Salpeter in form for any reasonable value 
for Merit (Pw0.19 for Merit=lOM0). Accordingly, we adjudge that 
there is not a deficit of white dwarfs, at least single objects, in this 
open cluster. If our conclusion is correct then these results suggest 
that the bulk of Praesepe white dwarfs did not gain a significant re- 
coil kick velocity from possible asymmetries in their loss of mass 
during the asymptotic giant branch phase of evolution. Fellhauer et 
al. (2003) find that if the mean kick velocity extended to a white 
dwarf during this phase is greater than twice the cluster velocity 
dispersion (the one-dimensional velocity dispersion of Praesepe is 
0.67±0.23kms~^ ; Madsen et al. 2002) then a significant fraction of 
these objects evaporate from the cluster (^60% at 600-700Myrs). 



4.2 The masses of the progenitor stars of the Praesepe white 
dwarfs 

The metalicities of the Praesepe and the Hyades open clusters 
are found, within uncertainties, to be very similar. For example, 
Cayrel de Strobel (1990) and Boesgaard & Budge (1988) determine 
[Fe/H]=0.10±0.06 and [Fe/H]=0.13±0.07 respectively for the for- 
mer, while Cayrel, Cayrel de Strobel & Campbell (1985) and Per- 
ryman et al. (1998) find [Fe/H]=0.12±0.03 and [Fe/H]=0.14±0.05 
respectively, for the latter. Furthermore, it has been recognised for 



decades that the space motions of the two clusters are comparable 
(e.g. Schwarzchild & Hertzsprung 1913, Eggen 1992, Madsen et 
al. 2002), leading to the suggestion that Praesepe is a member of a 
Hyades supercluster. If this is the case it is likely that the ages of 
the two clusters are similar. 

The bulk of determinations place the age of the Hyades cluster 
in the range 500-900Myrs (e.g. Barry et al. 1981, Kroupa 1995). 
Ferryman et al. (1998) have estimated an age of 625±50Myrs 
by fitting theoretical isochrones, which included a treatment for 
convective overshoot, to the Hipparcos-based cluster Hertzsprung- 
RusseU diagram. Claver et al. (2001) have compared modem the- 
oretical isochrones for ages 500, 700 and lOOOMyrs to the Prae- 
sepe sequence as observed in their V,V-I colour-magnitude dia- 
gram. They conclude that the upper main sequence is poorly repro- 
duced by the lOOOMyr isochrone but can be considered consistent 
with either of the two younger models. In addition, recent work on 
the X-ray properties of solar type Praesepe members reveals that 
these are similar to those of the F and G stars of the Hyades (Fran- 
ciosini et al. 2003). Furthermore, we note that the most detailed 
study to date of the spatial distribution and dynamics of the cluster 
(Adams et al. 2002) finds no evidence for the existence of a sub- 
cluster as suggested by Holland et al. (2001). Therefore, we follow 
Claver et al. (and more recently Ferrario et al. 2005) and in our 
subsequent discussion assume the age of Praesepe is 625±50M3Ts. 

We have derived the lifetime of the progenitor star of each of 
the four new white dwarf members by subtracting the estimated 
cooling time, shown in Table 1, from the adopted cluster age. The 
results of this process are shown in Table 2. For consistency, we ob- 
tained new spectroscopy of comparable quality for WD0837-I-218 
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and WD0837+185 (Figure 1) and analysed it using our new grid of 
model atmospheres. Thus we also provide in Tables 1 and 2, revised 
effective temperatures, surface gravities, masses, cooling times and 
progenitor lifetimes for these two objects. We note that the bulk 
of the discrepancy between the new and old parameter estimates 
for these two objects stems from the different datasets used in the 
two analyses. The spectral datasets in Dobbie et al. (2004) were 
of comparatively low quality and the true level of uncertainty in 
the associated effective temperature and surface gravity estimates 
appears to have been sUghtly imderestimated by our formal error 
analysis. 

Subsequently, we have used cubic splines to interpolate be- 
tween the lifetimes calculated for stars of solar composition by Gi- 
rardi et al. (2000) and have constrained the masses of these six pro- 
genitors to the values shown in Table 2. The quoted errors take into 
account the uncertainty in the cluster age. The locations in initial 
mass-final mass space of all six of these objects and the five orig- 
inal Praesepe white dwarfs, where their masses, cooling times and 
progenitor lifetimes and progenitor masses have been derived on 
the basis of the effective temperature and surface gravity measure- 
ments of Claver et al. (2001), are shown in Figure 3 (filled circles). 
We point out that five of the Praesepe white dwarfs virtually sit on 
top of each other in this plot. 

4.3 The initial mass-final mass relation 

Since the Hyades has a comparatively robust age determination we 
have added to Figure 3 the seven single white dwarf members of 
this cluster (open triangles). Their initial and final masses are de- 
rived from the effective temperature and surface gravity measure- 
ments listed in Claver et al. (2001). As there still appears to be 
considerable uncertainty as to the age of M37 (e.g. Kharchenko et 
al. 2005, Kalirai et al. 2005, Twarog et al. 1997), we do not include 
in Figure 3 recent data from this cluster. However, the eighteen 
white dwarfs from Praesepe and the Hyades only define the ini- 
tial mass-final mass relation for 2.7MQ^Mprogj54M0. Therefore, 
we have also added to this plot Sirius B (open circle; Liebert et al. 
2005a) and the white dwarf members of three other relatively well 
characterised but much younger open clusters (r^200Myrs), the 
Pleiades (open star), M35 (open diamonds; Williams et al. 2005) 
and NGC2516 (open '-l-'s; Koester & Reimers 1996). 

Using a new high S/N spectrum obtained with the WHT and 
our grid of TLUSTY models we have redetermined the effec- 
tive temperature and surface gravity of the only known Pleiades 
white dwarf WD0349+247 (LB1497) to be Teff=3284ltJ^^ and 
log g=8.63l jgg respectively (Figure 4). This corresponds to a mass 
of 1.02±0.02Mo, in excellent agreement with the gravitational 
redshift determination (l.O2±t^;^gM0; Wegner et al. 1991). For 
the white dwarf members of M35 and NGC2516 we adopt the ef- 
fective temperatures and surface gravities listed in Table 1 of Fer- 
rario et al. (2005). Excepting M35, where, given the result of Kali- 
rai et al. (2003), we prefer an age of 160±25Myrs, the progenitor 
mass for each of these white dwarfs is derived assuming the sys- 
tem/cluster age and metalicity adopted by Ferrario et al. (2005). 
The upper limits to the progenitor masses of the white dwarf mem- 
bers of NGC25 16 and the Pleiades displayed in Figure 3 have been 
determined by supplementing the model grid of Girardi with an 
evolutionary calculation for a 9M0 star drawn from a grid of a pre- 
vious generation of the Padova models (Bressan et al. 1993). 

An examination of the thirty stars in Figure 4 indicates that the 
majority appear to follow rather closely a monotonic relation be- 
tween their initial and final masses. Indeed, the white dwarf mem- 
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Figure 4. The results of our fitting of synthetic profiles (thin black lines) to 
the observed Balmer lines of WD0349+247, H-/3 to H-8, (thick grey Unes). 
The fluxx units are arbitrary. 



bers of the subsolar metalicity cluster M35 ([Fe/H]Ri-0.3; Sung & 
Bessell 1999) seem to form a natural extension of the relation de- 
fined by the Hyades and Praesepe white dwarfs (see Figure 3). This 
suggests that any effect metalicity may have on the form of the 
IFMR is probably not detectable at our current level of precision. 
Nevertheless, two stars, both attributable to Praesepe, appear to de- 
viate more significantly (<^3o-) from the general relation (LB5893 
and WD0837+218). As briefly discussed by Claver et al. (2001) 
and Dobbie et al. (2004) these white dwarfs appear to be too hot 
for their relatively high masses. We note here that their masses 
are coincident with the secondary peak in the mass distribution of 
white dwarfs (>O.8M0; see Figure 13 of Liebert et al. 2005b). As 
discussed by Liebert et al. it has been suggested that a significant 
proportion of the white dwarfs in this secondary peak may be the 
progeny of close binary systems. 

The initial mass-final mass data shown in Figure 3 can be 
reasonably approximated by a linear function. In deriving the co- 
efficients for this we exclude the two objects which appear to de- 
viate most strongly from the general relation as these most likely 
result from close binary evolution. Additionally, we choose to re- 
ject from the our fit the magnetic white dwarf EG61 as the true 
uncertainties in the effective temperature and surface gravity deter- 
minations for this star, obtained from the Zeeman split profiles of 
the Balmer lines, could be significantly larger than the formal er- 
rors quoted by Claver et al. (2001). Thus performing a linear least 
squares fit to the remaining twenty-seven white dwarfs we derive 
fit co-efficients ao=0.289±0.051, ai=0.133±0.015 (solid line). We 
note that had we adopted pure-C core white dwarf evolutionary 
models, the cooling times would have been systematically larger 
by between 5-30Myrs, resulting in shorter progenitor lifetimes, 
larger progenitor masses and a marginally flatter IMFR (ao=0.304, 
ai=0.127; dashed line). In either case our fit is somewhat steeper 
than the widely applied Weidemann (2000) relation, in particular at 
Mprog^4M0 where the latter flattens slightly (dotted line). 

In standard stellar evolutionary theory, the CO core of an in- 
termediate mass star can grow to ~1.1M0 before the carbon ig- 
nites. The resulting degenerate NeO core collapses, through elec- 
tron capture onto Ne, to form a neutron star. Extrapolating our fit to 
the initial mass-final mass data we find that a I.IM© white dwarf 
is produced by a star with an initial mass of ~6.1M0. However, 
more recent theoretical modelling by Garcia-Berro et al. (1997) 
have shown that it is possible for stable more massive NeO de- 
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generate cores to form via single star evolution. For example, their 
calculation of the evolution of a lOM© star ends in the production 
of a I.26M0 NeO white dwarf. However, a 1.37M© degenerate 
NeO core produced in a similar computation for a more massive 
star collapses, resulting in a type II supernova explosion. These 
findings led Weidemann (2000) to suggest that the upper limit for 
a white dwarf mass is ~1.3M0. On the basis of an extrapolation 
of our fit this corresponds to progenitor masses of 7.6(6.8-8.6)M0, 
somewhat lower than the 10-11 M© suggested by the models of 
Garcia-Berro et al. (1997). 

There is no observational evidence to indicate that the form of 
the IMFR can be represented by an extrapolation of our linear func- 
tion for Mprog^6M0 . Nevertheless, if we were to assume that such 
an extrapolation is valid and that the maximum white dwarf mass 
is I.3M0 then these results argue that the minimum mass of a core 
collapse supernova progenitor lies in the range '~6.8-8.6M0. We 
note that Ratnatunga & van den Bergh (1989) estimate a Galactic 
type II supernova rate of I.IIq g century"^ for a minimimi progen- 
itor mass of 8M0 . 



5 SUMMARY 

We have spectroscopically confirmed four more white dwarf mem- 
bers of Praesepe and argued that there are at least eleven white 
dwarfs in this cluster. We find that this number is consistent with 
that expected if the initial mass function of Praesepe was Salpeter 
in form. We suggest that this consistency indicates that few white 
dwarfs have evaporated from the cluster since their formation. If 
this is correct then any recoil velocity kick arising from possible 
asymmetries in mass loss during the final stages of stellar evo- 
lution is likely smaller than ~2.5kms~^. We find most stars ap- 
pear to follow, relatively closely, a monotonic relation between 
their initial and final masses. This relation can be reasonably ap- 
proximated by a linear function with coefficients ao=0.289±0.051, 
ai=0.133±0.015. Extrapolation of this linear function beyond ini- 
tial masses ~6M0 supports a minimum mass for a type II super- 
nova progenitor in the range ~6.8-8.6M0 for an assumed maxi- 
mtrni white dwarf mass of 1.3M0. 
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